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We consider a mixture of two species of spin-1 atoms with interspecies spin exchange, which may 
cooperate or compete with the intraspecies spin exchanges and thus dramatically affect the ground 
state. It represents a new class of bosonic gases differing from single-species spinor gases. We 
determine the exact ground states in several parameter regimes, and study the composite structures 
by using the generating function method generalized here to be applicable to a mixture of two species 
of spinor gases. The most interesting phase is the so-called entangled Bose-Einstein condensation 
(BEC), which is fragmented BEC with quantum entanglement between the two species, and with 
both interspecies and intraspecies singlet pairs. For comparison, we also apply the generating 
function method to a mixture of two species of pseudospin- 1 atoms, for which the total spin quantum 
number of each species is fixed as half of the atom number, in contrast with the case of spin-1, for 
which it is a variable determined by energetics. Consequently, singlet pairs in entangled BEC of a 
pseudospin-^ mixture are all interspecies. Interspecies spin exchange leads to novel features beyond 
those of spinor BEC of a single species of atoms as well as mixtures without interspecies spin 
exchange. 
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I. INTRODUCTION 

Spin-exchange scattering between bosonic atoms leads 
to novel ground states and interesting phenomena unex- 
plored previously [H-[E|. This remarkable subject broad- 
ens the scope of magnetism, which is traditionally based 
on spin exchanges of fermions instead. Mixtures of two 
species without a spin degree of freedom or, equiva- 
lently, mixtures of atoms of the same species with two 
spin states whose occupation numbers are both conserved 
have also been studied [6j . What about a mixture of two 
different species with interspecies spin exchange, in ad- 
dition to intraspecies spin exchanges? This question was 
first explored in a mixture of two species of pseudospin- 1 
atoms [7h1Q|| , where the most interesting phase was found 
to be BEC of interspecies singlet pairs, which was called 
entangled Bose-Einstein condensation (EBEC), or BEC 
with an entangled order parameter, emphasizing the as- 
pect that it is in an entangled state of two distinguishable 
atoms that BEC occurs. More generally, here we define 
EBEC as multispecies BEC with interspecies entangle- 
ment. It is a special kind of fragmented BEC. In EBEC, 
quantum entanglement, which is the most essential quan- 
tum feature not existing in classical physics, is amplified 
to a macroscopic quantum phase, just as BEC in a super- 
posed single-particle state amplifies single-particle super- 
position to a macroscopic phase, leading to the Josephson 
effect. EBEC also bears some similarities to the low- 
est energy state of a SU(2) symmetric model of a single 
species of pseudospin- ^ atoms which are forced to occupy 
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two orbital modes by conservation of the total spin in the 
cooling process [ll|, but there are also differences [lfjj . 

To pick out two hyperfine states of an alkali atom 
representing pseudospin- i avoiding spin-exchange loss 
to other hyperfine states may require some careful tun- 
ing In addition, with the degree of freedom being 
pseudospin, there is no reason to expect rotation sym- 
metry in the interaction. In contrast, in a mixture of 
two species of spin-1 atoms with full access to the spin- 
1 multiplet, the total single-particle energy of any two 
scattered particles is always conserved, whether in the 
absence or the presence of a magnetic field, hence spin- 
exchange scattering is energetically protected. More- 
over, the total spin of each species is conserved by 
the intraspecies interaction, while the total spin of the 
whole mixture is conserved by the interspecies interac- 
tion. Therefore, if it is also a ground state of a mix- 
ture of two species of spin-1 atoms, EBEC may be more 
experimentally accessible and stable in such a mixture 
than in a pseudospin- ^ mixture. However, our consid- 
eration of spin-1 mixture had been impeded by the ap- 
parent complexity of the spin-exchange interaction be- 
tween two spin-1 atoms of different species, until it was 
shown recently that it is simply of Heisenberg form [ljj ■ 
Some mean-field-like investigations of spin-1 mixtures 
have been carried out, but possible entanglement between 
the two species was ignored [ill EH ■ 

In this article, we rigorously study the ground states 
of a mixture of two species of spin-1 atoms in various pa- 
rameter regimes. Bosonic symmetry within each species 
and its absence between different species together lead 
to rich structures, with interesting features beyond those 
of a single species of spinor gas. For a spin-1 mixture, 
we find EBEC in certain parameter regimes, with the two 
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species significantly entangled. Two atoms of the two dif- 
ferent species can form interspecies singlet pairs, with dif- 
ferences from a pseudospin-^ mixture, however; in EBEC 
of a spin-1 mixture, intraspecies and interspecies singlet 
pairs coexist. 

The rest of the paper is organized as follows. The 
many-body Hamiltonian is given in Sec. II. Under a com- 
mon assumption for spin-1 bose gases, the present Hamil- 
tonian can be written solely in terms of spin operators of 
the two species and of the total system. Then in Sec. Ill, 
we find the ground states in terms of spin quantum num- 
bers, in various parameter regimes. These ground states 
are given in terms of boson creation operators in Sec. IV. 
In Sec. V, composite structures of these ground states 
are studied by generalizing a generating function method 
from a single species to a mixture. In Sec. VI, we discuss 
a mixture of two species of pseudospin-i bosons. The 
paper is summarized in Sec. VII. 

II. THE MANY-BODY HAMILTONIAN 

For two spin-/ atoms of different species, there is no 
permutation symmetry between them, hence the total 
spin can be F = 0, • • • 2/. The effective interaction is 
thus 

2/ 

V(r a - r b ) = 5(r a - r b ) £ gfP F (1) 

F=0 
2/ 

= *(r -r 6 )5^(F .F 6 )', (2) 
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where g F h is the interaction strength proportional to the 
F-channel scattering length, and Pp is the projection 
operator for the total spin F, and can be expanded in 
terms of 1, F a • F b , • • • , (F„ • F b ) 2 '. For / = 1, cf = 
-g- b /3 + gf + gf/3, cf = -gf/2 + gf/2, and cf = 
5o7 3 - 9i V 2 + ff/6- It has been shown that cf = 
3gr/4 + gs/4, cf = (g T - ffs)/16, and cf = 0, where g T 
and gs correspond to the triplet and singlet states of the 
two valence electrons of the scattering atoms [12j . 
Therefore the many-body Hamiltonian is 

H = H a + n b + H ab , (3) 

where 

+ i f driplt^pffiS^Spo- + c%F a ^ ■ F QpCT )Vw!/w 

(4) 

is the well-known Hamiltonian of spin-1 atoms |l| of 
species a (a = a,b), while 

Hab = J dr^ M V4 P ( 5 o b( W < W + cfFa^ ■ F bp(T )tp b<T tp au 

(5) 



is the interspecies interaction. Here the field operator 
ipaf* corresponds to spin /i component of species a (fx — 
— 1,0,1), and F QMI , represents the (/xz/) element of the 
spin-1 matrix of species a. 

h 2 

h a = - — V 2 + U a {r)- la B-F a (6) 

is the single particle Hamiltonian of species a, m a and 
7 a are the mass and the gyromagnetic ratio of an atom 
of species a, respectively, B is a uniform magnetic field, 
and U a is the trapping potential for an atom of species 
a. 

For the single-particle orbital wave function, certainly 
one can use the usual single-mode approximation to write 
"0QM( r ) = a M^( r )> where a M = a M ,fo M is the annihila- 
tion operator and (f> atl is the lowest single-particle orbital 
wave function for species a and spin /i. Then we obtain 
a Hamiltonian in terms of creation and annihilation op- 
erators, with integrals of various products of <^ QjU 's and 
their complex conjugates entering as coefficients of the 
terms of the Hamiltonian. 

Nevertheless, to simplify the matter, we can follow 
an additional common assumption for spin-1 bose gases, 
namely that the single particle orbital wave function is 
mainly determined by the spin-independent part of the 
Hamiltonian and is thus independent of spin; that is, 
<fian = 0q is independent of spin /i. For a homogeneous 
system, <p a = l/VO, where 51 is the volume. But our 
discussions also apply to the inhomogeneous case. 

Therefore, the Hamiltonian can be simplified as 

where 

S Q = a^F^au (8) 

is the total spin operator for species a, "f a ,Jb > 0, 
C = N a e a + N b e b + f {N 2 a - N a ) + C ^{N 2 - N b ) - 
c a N a -c b N b + cfN a N b is a constant, c£ = eg / d 3 r|0 Q | 4 , 
c t b = I d 3 r\4> a 4> b \ 2 , (k — 0,2), and we have simplified 
notations and cf as c a and c ab , respectively. 

If c ab = 0, then there is no entanglement between the 
two species of atoms, the ground state is simply a direct 
product of the ground states of the two species of spin-1 
atoms. 

We now set out to find the ground states of ([7|) in 
various parameter regimes. 



III. GROUND STATES IN TERMS OF SPINS 

We assume 7 a = 7^ = 7 > 0, as satisfied by alkali 
atoms with the same nuclear spin, for example, 7 Li, 23 Na, 
39 K, 41 K and 87 Rb, all of which have nuclear spin 3/2. 
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Then S a and Sb together with the total spin S and its 
z-component S z are all conserved. Therefore, the ground 
state is 

\G) = \S a ,S b ,S,S z ), (9) 

with the four spin quantum numbers being the inte- 
gers that minimize the energy E = c ~ c S a {S a + 1) + 

£=f^S b (S b + 1) + ^S(S + 1) - -fBS z , where B > 
is the magnitude of the magnetic field and the constant 



In the following, we find the ground states in the 
form of (|10|) . by minimizing for various parameter 
regimes. For specification, we assume N a > N b without 
loss of generality. 



A. c ab < 

First, we consider the cases with c ab < 0. Then for 
given S a and Sb, it is S = S z = S a + S b that minimizes 
the energy (fTT|) . 

The ground state is thus 

\G(c ab < 0)) = \S™,Sl n ,S™ +S?,S? +SD, (12) 

where <S™ and S™ are respectively the values of S a and 
S b that minimize E{c ab < 0) = ^f^S a {S a + 1) + 

S^ls b (Sb + l) + ^(S a + Sb)(S a + S b + l)-~fB(S a + Sb), 
and must be calculated separately for different subcases. 

In particular, here we consider c ab < 0, c a < c ab and 
c b < c ab . E always decreases as S a or S b increases. For 

I 



where S z = —N a — Nb, ■ ■ ■ , N a + Nb, the prime in the 
subscript to \G)p represents the absence of a magnetic 
field, and gi'(Sbz) is the Clebsch-Gordan coefficient. 

With c ab < 0, one can see that the ground state re- 
mains as |G)j or \G) r if c a = c ab while c b < c ab , or 

c b _ c ab while c a < c ab^ QJ . c a _ c b _ c ab Thig ig be _ 

cause to minimize E, first S is maximized to the largest 
possible value N a + N b . Consequently S a and S b have to 



C has been neglected. Note that the minimization of E 
is under the constraints \S a — Sb\ < S < S a + S b and 
-S <S Z < S. 

Moreover, for a given S, S z = S minimizes the energy. 
Therefore the ground state must be 

\G) = \S™,SF,S m ,S m ), (10) 

where 5™, <S™ and S m are, respectively, the values of S a , 
S b and S that minimize 



(11) 

I 

each species of spin-1 atoms, the largest possible value of 
the spin quantum number is clearly its number of par- 
ticle. Hence 5™ = N a and 5 b m = N b . Therefore the 
ground state is 

|G>j = |JV a ,JV6,JV + JV6,JV a + JV 6 ) (13) 
= \N a ,N a ) a ®\N b ,N b )b, (14) 

where the subscript "I" in \G)i denotes the parameter 
regime. Similar notations are used for ground states in 
other parameter regimes. Throughout this paper, a spin 
basis state written without a subscript is of the total 
system, while a spin basis state with a subscript "a" is 
of species a; that is, 

I S a , S az ) a , 

is a state of species a with total spin S a and its z com- 
ponent S z . Obviously \G)i is the direct product of the 
ferromagnetic states of the two species. 

In the case of c ab < 0, c a < c ab and c b < c ab , if the 
magnetic field is absent, that is, B = 0, then there are 
2(N a + Nb) + 1 degenerate ground states 

(15) 
(16) 

I 

be maximized to N a and Nb, respectively, even though 
one or both disappear in E. Therefore regime I can be 
expanded as c ab < 0, c a < c ab and c b < c ab . 



I 

„a ab „b ab ab 

E = S a (S a + 1) + — ^— S b (S b + 1) + —S(S + 1) - 1 BS. 



\G) r - \N a ,N b ,N a + N b ,S z ) 

= 9i>(S bz )\N a ,S z -S bz ) a ®\N b ,S bz ) b , 

S bz = -N b 
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B. c ab > 

Now we turn to cases with c ab > 0, for which it is 
useful to rewrite E as E(c ab > 0) = ^f^S a {S a + 1) + 
^S b (S b + l) + ^(S+±-^r-^-Z§r. When 
c ab > 2yB, c a - c ab > 0, and c b - c ab > 0, the ground 
state is 



\G) n = |0,0,0,0) (17) 
= |0,0) o ®|0,0> 6) (18) 



|G) //J = |JV,iV,0,0) 
1 N 

= 1 y 

2iV + 1 ^ 

m— — 



which is a maximally entangled state, as one can easily 
see by considering the reduced density matrix of either 
species. 



which is a maximally entangled state. 

If c ab > 2-fB > or c ab > 2jB = under the con- 
straint N a = N b = N, one can see that the ground state 
remains as \G)m if c a < c ab while c b = c ab , or c 6 < c ab 
while c a — c ab . This is because to minimize E, S is 
minimized to while S a or Sb is maximized to N and, 
consequently, the other Sb or S a is maximized to N also, 
although it does not appear in E. 

If c ab > 2yB > or c ab > 2jB = 0, while c a = 



where giv(Sbz) is the Clebsch-Gordon coefficient. 

If < c ab < 2-/B, and c a - c ab < while c b = c ab , then 
the ground state is \N a , Sb,n,n), under the constraint 
l-^o-^l < n < N a + S b . If < c ab < 2 1 B, and c b -c ab < 



which is the direct product of the two singlet states of the 
two species. The result is valid even when c ab = 2-yB = 0, 
as S a and S b are still minimized, respectively, although 
such a trivial case is not our focus. 

If c ab > 2jB > or c ab > 2-fB = 0, one can see 
that the ground state remains as \G)u if c a > c ab while 

c b = c ab^ Qr c b > c ab wmle c a = Q ab _ Tms j g becauge 

to minimize E, S and S a or S b are minimized to and, 
consequently, S b or S a has to be too, although it does 
not appear in E. 

Now let us consider the case with c ab > 2jB > or 
c ab > 2-fB = 0, together with conditions c a - c ab < 0, 

c b _ c ab < q and Na = Nb = N _ ground gtate ig then 

the global singlet state 

(19) 

l) m \N,m) a ®\N,-m) b , (20) 

I 

In general, for two species of spin-/ atoms of an arbi- 
trary /, with N a = N b = N, a singlet state with maximal 
S a and S b is 



(21) 



I 

c b = c ab , then the ground state is degenerate, and is 
156,56,0,0), with < S b < N b . 

Finally, if < c ab < 2jB, c a -c ab < 0, and c b -c ab < 0, 
while \N a - N b \ < n < N a + N b , where n = Int($ - \), 
then E is minimized when S a — N a , Sb — Nb, and S = n. 
Here Int{x) represents the integer closest to x and in the 
legitimate range given above. is the integer no larger than 
and closest to x. Hence the ground state is 

(22) 
(23) 

I 

while c a = c ab , then the ground state is \S a ,N b ,n,n), 
under the constraint \S a — N b \ < n < S a + N b . If < 
c ab < 2-yB while c a = c b = c ab , then the ground state 
is \S a ,S b ,n,n), under the constraint — S b | < n < 



J 



1 fN 

\fNJN,0,0)= y {-l) m \fN,m) a ®\fN,-m) bl 

•' ' m=-/JV 



I 

|G>/v = \N a ,N b ,n,n) 

= y 9iv(S bz )\N a ,n- S bz ) a ®\N b ,S bz )b, 

S b z = -N b 
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Sa + S b . 

Note that the preceding ground states are all unique, 
because in each case, not only S and S z , but also S a 
and S b are specified. It has been known that for a single 
species of spin-1 atoms, \S a ,S a } a is unique [14|. Hence 
any spin basis state 



\S a ,S az ) a = {r(S a )---r(S az +l)}- 1 (S a _) s * 
is also unique, where 

S a _ = V2(a\ao + aja_i) 
is the spin lowering operator, 

r{m) = y/(S + m)(S -m + 1). 



| S a ) $ a ) O 

(24) 
(25) 

(26) 



species of spin-1 atoms, as can be obtained from the spin 
basis states of the two species, is unique. 



IV. GROUND STATES IN TERMS OF BOSON 
OPERATORS 



We now proceed to determine the expressions of these 
ground states in terms of boson operators. First, it is 
straightforward to obtain 



\G)i = \N a , N b , N a + N b , N a + N b ) 
= , \ , (4) Na (b\) Nb \0), 



Therefore any state \S a , S bl S, S z ) of a mixture of two and 
I 



(27) 
(28) 



\G)n = |0, 0,0,0) 

= Z u [2a\aU (al)Y° /2 [2b\bU - (bl)Y b/2 \0), 



(29) 
(30) 



where Z u = [(N a /2)l(N b /2)l2^+ N ^/ 2 (N a + l)\l(N b + 
l)!!] 1 / 2 is the normalization constant Q. From this ex- 
pression, it can be seen that \G)j exists only if N a and 
N b are both even. 

One can rewrite in terms of boson operators \G)p, 
\G)in and |G)iy, in which cases S a takes the largest 
possible value N a , by calculating the spin basis states of 
each species using (|24"1) . with 



I AT N ) 



|0). 



(31) 



I 

and then substituting them into (jTHJ) , (|2"U| and (f2"3"|) . re- 
spectively. For convenience in reading, we write them 
down explicitly in the following equations. 



\G)i> = \N a ,N b ,N a + N b ,S z ) 

N b 

- E 

S bz =-N t 



9v{S bz )2 ^~ 



^NjNj-r(N a ) ■ ■ ■ r(S z - S bz + l)r(N b ) ■ ■ ■ r(S bz + 1) 



x (ajai + aLxao) 



Na-S.+S b ,J Na flf 



'aribl^ + bUbar^bf'lO), 



(32) 



(33) 



where 



gr(S bz ) 



' (2iV a )!(2jV b )!(iV a + N b + S z )\(N a + N b - S z )\ 

(2N a + 2N b )l(N a -S z + S bz )l(N a + S z - S bz )l(N b - S bz )l(N b + S bz )V 



(34) 



\G)m = \N,N,Q,0) 
2 N 



(2N + 1)N\ 



N 

E 



(-iy 



-N 



r(N) ■ ■ ■ r(m + l)r(N) ■ ■ ■ r(-m + 1) 



x(4ax +aL 1 ao) JV - m « t 1 Ar (^i +bl 1 b o ) N+m bf\0), 



(35) 



(36) 
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\G}iv = \N a ,N b ,n,n) 



N a + N b -7, 



E 



giv(S bz )2 2 



S bz = -Ni 



JNjNj.r(N a ) ■ ■ ■ r(n - S bz + l)r(N b ) ■ ■ ■ r(S bz + 1) 



x (ajai +- a^ao) 7 



10), 



where 



giv(S bz ) = (-1) 



' (JVg + AT b - n)!(2n + l)!(iV a + n- S bz )l(N b + S bz )\ 

(N a +N b +n+ l)!(n + AT a - - ^ a + JV 6 )!(iV a _ „ + S te )!(jV 6 - 5 bz )! " 

I 



(37) 



(38) 



(39) 



1 5;,, 5(,, 0, 0), which is the ground state in some bound- 
ary regime of regime III, can be given by (|36[) . with N 
replaced by S b . \S a , S b ,n,n), which is the ground state 
in some boundary regime of regime IV, can be given by 
([55)) and dSH]), with A" a and N b replaced by S a and S b , 
respectively. 



A clearer picture of the composite structures of these 
ground states in terms of basic units is revealed by gen- 
eralizing the method of generating function 0, [l4[ , in the 
next section. 

I 



V. COMPOSITE STRUCTURES OF THE 
GROUND STATES 

A. Generating function method for \S a , Sb, S, S) 

Now we consider the construction of |5 ,5ft, 5, 5) of 
a mixture of two species of spin-/ atoms, with total 
spin and its z-component both being 5. With S z being 
maximal, we may consider a configuration of the state 
|5 a ,5ft,5, 5), in which there are Q mj ,n<,l< copies of unit 
j, which is made up of rrij a-atoms and rij &-atoms and 
carrying spin lj . Denoting the creation operators for unit 

3 as @L. n.-.L) we have 



\s a ,s b ,s,s) = ^2A({Q mjinjilj })H(el jtnjtlj ) 

I 



Q.,..i. 



10), 



(40) 



where A is a coefficient and the summation is over all 
possible values of {(),„ „ ./ | . 

The possible values of //lj, Tij lj and Qnij ,rij ,ij are 
subject to the constraints from N a , N b and 5, 



find. 



H ,(}:, ,n ./ = N a , 

3 

lllj(.>,:. ,« ./ = N b) 

J 

y, tjQ'" " / = s> 



(41) 



For an integer /, we define the generating function as 



G(x a ,x b ,y) = ]T M{N a ,N b ,S)x^x^y s , (43) 

N a ,N b ,S 

where x a , x b and y are complex numbers inside the unit 
circle and M(N a , N b , 5) is the number of solutions of the 
sets of the nonnegative integers {Qm,j,n 3 -,lj}- Following 
the method of [14]], we obtain that 



as well as the constraints from S a and S b , 

5^|5 a , 5;,, 5, 5) = 5 a (5 a + l)|5 a , 5ft, 5, 5), 
Si \S a , 5;,, 5, 5) = 5fc(5fc + l)|5 a , S b , 5, 5). 



(42) 



where 5^ = {a\a\ — ali^-i) 2 + 2(aJ) 2 aia_i +- 
2a 1 a t _ 1 aQ + 2aQaoaia\+2al l aoa-ia ji _ 1 , as one can easily 
I 



G(x a ,x b ,y) = 
dz I — z^ 1 



n n a 



c 2ni z-y 11 (l-x a z^)(l-x b z^)' 

3a — J 3b — / 



where the contour integral is along the unit circle C. 
For / = 1, we obtain that 



(44) 



G{x a ,x b ,y) = J2( x a 



i,i,o+2Q2,o,o+Q i,o,i Qi,i,o+2Q ,2,o+Qo,i,i ..Q^o.i+Qo.i.i 



ft V 

E Qi.i,o+2Q2,o,o+Qi, o.i + l Qi,i.q+2Qo.2,o+Qq,i,i+1 Q 1 1 +Q i i + l \ 
x « x ft y )■ 



(45) 



where the summations are over all possible values of nonnegative. Comparing (|43[) and (|45[) . we have 



Qi,i,o, Q2,o,o, <9o,2,o, Qi,o,i and Qo,i,i, all of which are 



M(N a , N b , 5) = Mi(N a , N b , 5) + M 2 {N a , N b , 5), (46) 
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where Mi(N a , N b , S) is the number of solutions to the 
set of equations 



Qi,i,o + 2Q2 : o,o 

Ql,l,0 + 2Qo,2,0 

Qi,o,i 



Ql,0,l = A a , 

Qo,i,i = A&, 

Qo,l,l = "5; 



(47) 



while M2(N a , Nb, S) is the number of solutions to the set 
of equations 



(48) 



Ql,l,0 + 2<32,0,0 + Ql,0,l + 1 = N a , 

Qi,x,o + 2Qo,2,o + Qo.i.i + 1 = Nb, 
Qi.o.i + Qo.i.i + 1 = S. 



In general, there may be multiple solutions to (|47p . In 
each solution, there are Qi,i,o interspecies singlets con- 
sisting of one a-atom and one 6-atom, Q2.o,o singlets con- 
sisting of two a-atoms and Qo,2.o singlets consisting of 
two o-atoms. In addition, there are Qi,o.i a-atoms with 
z-component spin 1, as well as Qo,i,i fr-atoms with z- 
component spin 1. 

There may also be multiple solutions to (14"51) . In each 
solution, there are Qi,i,a interspecies singlets consist- 
ing of one a-atom and one 6-atom, Q2,o,o singlets con- 
sisting of two a-atoms and Qo,2.o singlets consisting of 

I 



two 6-atoms. Also, 
z-component spin 
component spin 0, 
z-component spin 
z-component spin 

component spin 0, as well as Qi,o.i 
z-component spin 1. 



either there are Qi,o,i 
1 together with one a- 
as well as Qo.i.i + 1 
1; or there are Qo,i,i 
1 together with one b- 
1 



a-atoms with 
atom with z- 
6-atoms with 
6-atoms with 
atom with z- 
a-atoms with 



B. 



\G)i and \G)u 



As a simple example, we can verify that |G)j = 
\N a ,N b ,N a + N b ) = \N a ,N a ) a ® \N b ,N b ) b is indeed as 
given in (j2"5j) . For S = N a + N b , the only solution to 
63 is Qi,o,i = N a , Q ,i,i = N b , and Qi,i,o = Q 2 ,o,o = 
Qo.2,0 = 0, while there is no solution to (|4"5|) . Hence in 
this case, \G)i must be given by (j2"5j) . which is obviously 
an eigenstate of S% and S b , as it should be. 



For S = 0, there is no solution to (|4*5|) , while there 
are solutions to (|47|) with Qi,o,i = <3o,i,i = 0, Q2,o,o = 
(N a - Qi,i,o)/2, and Q Q . 2 ,o = (N b - Qi,i,o)/2- Hence a 
state \S a , S b , 0, 0) can be expressed as 



-Qi,i,o)/2/ot 



(e5,2,o) 



,o)/2 



10), 



(49) 



where 



t _ 



= a\b t _ 1 - 



,t 6 t 



-A 



(50) 



is the creation operator for an interspecies singlet pair, 

while B 2 Q q = 2a}a_ 1 — a and O 2 = 2b[b_ 1 — b are 
creation operators of intraspecies singlet pairs of a-atoms 
and o-atoms, respectively. \S a , S b , 0, 0) is a superposition 
of configurations with all possible values of Qi,i,o ; i n each 
of which there are Qi,i,o interspecies singlet pairs, (A a — 
Qi,i,o)/2 singlet pairs of two a-atoms and (N b — Qi,i,o)/2 
singlet pairs of two 6-atoms. In addition, (j4"2")l with S a — 
S b = N imposes constraints determining C(Qi,i.o)- 

Equation {HHJ indicates that in \G) H = |0, 0) a (g)|0, 0) b , 
the two species are disentangled. Therefore Qi,i,o = 0, 
and thus Q 2 ,o,o = N a /2, and Q ,2,o = N b /2; that is, 
there are only intraspecies singlet pairs. Hence ([50)1 is 
confirmed. As mentioned previously, this is subject to 
the condition that N a and N b are both even. 

Similarly, we would like to note that for a gas of a 
single species a of A spin-1 atoms with c a > 2-fB, if N is 

I 



I 

odd, the (non- normalized) ground state should be |1, 1) = 

al(2alaLi-c4 2 ) (Ar ~ 1)/2 |0), rather than the singlet |0,0), 
which exists only for even A. 

C. \G)m and \G) IV 

Now we look at \G) IH = \N, A, 0,0). Equation (dUJ) 
already indicates that the two species are strongly entan- 
gled, therefore \G)ni is in the form of (|49l) . 

For \G)iv, one considers d47j) and (gSJ with S — S z — 
n. It can be found that the solution to (j4"T|) is 



N a -N b +n 



Qi,o,i 
Qo.i,i 



which is valid if A a + N b 
to HI is 



Q2,0,0 
<3o,2,0 



(51) 



^2,0,0, 

n is even, while the solution 



Qi,o,i 

Qo,i,i 



N a -N b +n-l 



N, 



-iV?+n-l 



Q2,0,0 — Qo,2,0, 
Qo,2,0 ~ Q2,0,0- 



(52) 



which is valid if A a -f 
Hence, if A a + A 6 



N b — n is odd. 
- n is even, 



\G)iv = \N a ,N h n,n) 



^Wl,l,0iV2,0,0,V0 ; 2,0j fl i ° 1 "1,1,0 
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2,0,0 



0,2,0 



''°' 2 >>, 



(53) 
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where Qi.0,1 and Qo.1,1 are given by (|5"Tj) while if N a + N b — n is odd, 



\G)iv = \N a ,N b ,n,n) = J2 A (Qi,i,o,Q2,o,o,Qo,2,o)ala\ Ql} °' 1 b\ Q ° }1 ' 1 + @I,i, Qll ' ©2,o,o w o._>j> 

1„| Q2,0,0„| C?0,2,0 

"2.0.0 "O^.O 



■Ei4 / (Qi,i,o,Oa,o,o.Qo,a,o)oi &5 & i e i,i,o 



|o>, 



(54) 



where Qi,o,i and Qo,i,i are given by ([52]) . In both (|53|) 
and (I54|) . the summations are over Qi,i,o> Q2,o,o and 
Qo.2,0- The coefficients A and A' are determined by con- 
straints gl) with S a = N a and S b = N b . 



VI. COMPOSITE STRUCTURES OF |5, S> OF A 
PSEUDOSPIN-1/2 MIXTURE 

As a comparison, we now apply the generating function 
method to determine \S,S) of a mixture of two species 
of pseudospin-1/2 atoms. Unlike the case of spin-1, for 
each species a of pseudospin-i, the total spin is always 
fixed to be S a — N a /2. With N a and N b fixed, the spin 
state of the pseudospin-^ mixture is only determined by 
5 and S z of the total spin; that is, 



15,5, 



\\N a ,±N b ,S,S, 



(55) 



The basic method for constructing the maximally po- 
larized 15, S) remains the same as that for a spin-1 mix- 
ture, as described in Sec. V.A. In a configuration of the 
state 15 5), there are Q m ,n .1 copies of unit j, which is 
made up of rrij a-atoms and rij 6-atoms and carries spin 
lj/2. The generating function is now defined as 

G(x a ,x b ,y) = J2 M(N a ,N b ,S)x^x b Nb y 2S . (56) 

N a ,N b ,S , 



I 

It can be obtained that 

G(x a ,x b ,y) = 



[ ^ l ~ z TT TT - 

J c 2m z-y 11 f 11 f (1 - x a z 2 ^)(l - x b z 2 ^)' 



-2 / / 

3a=~f 3b=~f 



where the contour integral is along the unit circle C . Then 
one obtains 



G(x a ,x b ,y) 
Therefore, 



_ \ -Oi.i.o+Qi.o.i-Qi.i.o+Qo.i.i, Qi, ,i+Qo,i,i 



v 



Qi,i,o 
Qi,i,o 
Qi,o,i 



Qi,o,i 
Qo,i,i 



N a , 

N b , 
25 



(58) 



(59) 



which has a unique solution, 

Qi,o,i = S+UN a ~N b ), 

Qo.i.i = S+hN b -N a ), 

Qi.i.o = S+l(N a + N b ). 
which is subject to the condition that Qi,i,o 5 Qi,o,i and 
Q '0,1,1 should all be nonnegative integers. Under this con- 
dition, 



(60) 



|5 5) = a. 



r S+(N a -N„)/2 ,S+(N b -N a )/2 



$1,1,0 



Jj S+(N a +N b )/2 



10), 



(61) 



where the normalization constant is neglected, and 

$ 1)1]0 t = a \b\ - a\b\ (62) 



is the interspecies singlet creation operator for a 
pseudospin-i mixture. A sufficient and necessary con- 
dition for Qi^i^, Qi,o,i and Qo,i,i all to be integers is 
that 25 + N a + N b is an even integer. 

When 5 = 0, the only consistent solution of is 

I 



r 



Qi,o,i = Qo,i,i = while Qi,i,o = N if and only if N a 
N b = N. The state is then the global singlet 



|0,0) = 0>i,i,o tAr |0) 



(63) 



When N a > N b , the smallest value of 5 is (N a - N b )/2, 
for which Qi^o = N b , Qi, ,i =N a - N b , and Q ,i,i = 0, 
hence 



N a -N h 



\±(N a - N b ), -- (N a - N b )) = a\" a " b ^ Nb \0) 



(64) 



of which the singlet state (jSBf is a special case. 



VII. SUMMARY 



To summarize, we have considered a mixture of two 
different species of spin-1 gases with interspecies spin- 
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No. 


Parameter 
regimes 


Ground states 


T 
1 


c ab < 0, 

„a , ab 
C < C , 

c b < c ab . 


\G)i = \N a , N b , N a + N b , N a + N b ) 

1 AT AT \ K> 1 AT AT \ 
= \J\a,Na)a <8> | JVb, iVj,)j, 

(disentangled) 


II 


c ab > 2jB, 
c a > c ab , 
c b > c ab . 


|G)zi = |0, 0,0,0) 
= |0,0)„ ® |0,0);, 
(disentangled) 


III 


c ab > 2^B > 
or c ab > 2-yB = 0, 

a ab 
C < C , 

C < C 

N a = N b = AT. 


|G)/// = |JV, JV, 0, 0) 
(entangled) 


IV 


< c ab < 27B, 
c a < c a6 , 
< c ai \ 


\G)iv = \N a ,N b ,n,n) 
n = Int{^-\) 
(entangled) 



TABLE I. Ground stats of a mixture of two spin-1 atomic 
gases in four parameter regimes. The first (I) is a direct 
product of ferromagnetic states of independent BEC of the 
two species. The second (II) is a direct product of the two in- 
dependent singlet states of the two species. The third (III) is 
an interspecies singlet state, which we call EBEC. The fourth 
(IV) is also EBEC if n < N a + N b . Exactly speaking, \G)n 
is subject to the condition that N a and N b are both even. 
Situations on the boundaries of these parameter regimes are 
discussed in the text. 



exchange scattering, as an extension of our previous work 
on a mixture of two different species of pseudospin- i 
gases, going beyond the usual spinor bose gases and 
BEC mixtures without interspecies entanglement. Inter- 
species spin exchange favors spin ordering between differ- 
ent species, while intraspecies spin exchange favors spin 
ordering within each species. The ground state of such 
a mixture thus depends on the parameters in the many- 
body Hamiltonian, which we have shown to be reduced 
to a Hamiltonian of two giant spins. 

We have worked out the ground states in four typical 
parameter regimes, which are now reported in Table I. It 
is straightforward to verify that they are all fragmented 
BEC, by calculating one-particle reduced density matri- 
ces. When c° and c b are both less than c ab , which is 
negative or 0, all atoms of each species form a ferromag- 
netic state with the spin of each atom being ji = 1; that 
is, the ground state of the mixture is the direct prod- 
uct of two independent ferromagnetic states. When c a 
and c b are both larger than c ab , which is larger than 
or equal to 2jB, the atoms of each species form a sin- 
glet state, with the total spin of each species being 0; 
that is, the ground state of the mixture is the direct 
product of two independent singlet states, subject to the 
condition that N a and N b are even (otherwise there is 
minute deviation). These two ground states are disen- 
tangled between the two species. For N a — Nt = N, 
when c ab > 2^B > or c ab > 2jB = 0, while c a and 
c b are both less than c ab , the ground state is a global 
singlet state \G)ni = |iV, iV, 0, 0) , with total spin zero. 



When < c ab < 27_B,and c a and c b are both less than 
c ab < 0, the ground state is \G)iv = \N a , N b , n, n), where 
„ = /nt(2§-!). 

The latter two ground states exhibit EBEC, display- 
ing strong interspecies entanglement. A consequence of 
this entanglement is that the particle number in each 
spin state of each species is subject to strong quantum 
fluctuation. There are rich composite structures due to 
interspecies entanglement. By using the generating func- 
tion method, it has been revealed that \G)m and |G)jy 
are each superpositions of configurations with both in- 
traspecies and interspecies singlet pairs. It is interesting 
to note that \G)m = \G) IV when c ab = 27B, N a = N b , 
implying that \G)m and |G)/y belong to the same quan- 
tum phase. 

We have also used the generating function method to 
find the spin state \S, S) of a mixture of two species of 
pseudospin-i atoms with interspecies spin exchange. As 
the total spin of each species of pseudospin- \ atoms is 
always half of the atom number, the composite structure 
of a pseudospin- i atoms is simpler, with only one config- 
uration. Consequently, there are only interspecies singlet 
pairs when the total spin of the mixture is zero, hence 
EBEC in such a case is simply BEC occurring in an inter- 
species singlet state. Such a simplicity is lost in a spin-1 
mixture, in which intraspecies singlet pairs coexist with 
interspecies singlet pairs, and EBEC is generally defined 
as BEC with interspecies entanglement. Previous stud- 
ies on pseudospin- i mixture demonstrated that EBEC 
leads to various physical properties different from those 
of usual BEC, to be similarly studied in spin-1 mixtures. 
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Note added: Recently we became aware of a paper [l6[ 
that simply gives \N a , Nb, N a — N b , S z ) as the ground 
state "for large antiferromagnetic spin-exchange interac- 
tion between the two species" , without actually show- 
ing it gives the lowest energy. We had determined that 
the ground state varies with the parameters, and is not 
always this state even for large antiferromagnetic inter- 
species spin-exchange interaction. They also applied the 
same generating function method as ours to the Hamil- 
tonian with an additional Pq term, which, to our under- 
standing, had been shown to vanish by Luo et al. [l2j]. 
We also disagree with their statement that there is no 
interspecies singlet pairing when this additional Po term 
vanishes. We have worked out all possible ground states, 
which will be discussed elsewhere 
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